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1. Introduction 

■ The evaluation of the thermal average of physical quantities is one of the 

, main aims in statistical mechanics. The density matrix of a system is the 

most fundamental quantity to achieve this aim. Its diagonalization, how- 
ever, becomes exponentially difficult with growing system size L. One in- 
evitably has to give up this procedure in the thermodynamic limit. An 
alternative approach for quantum systems is to diagonalize the Hamilto- 
nian, and to sum up the contributions from each eigenstate. This means to 
divide the problem into two parts: (1) diagonalize, and (2) sum up. Again, 
both procedures become exponentially difficult with the increase of L. 

In this article we re-consider this problem for integrable quantum spin 
chains. We will show how the integrability helps bypassing the difficulties 
and yields exact estimates. The first problem, the diagonalization of the 
Hamiltonian, can, in principle, be solved by the celebrated Bethe ansatz. 
The second step, however, remains as a cliff wall. A first breakthrough, the 
string hypothesis approach, was achieved in the early 70's. 1,2 In this ap- 
proach one introduces so-called root density functions of strings and holes 
of various lengths for the diagonalization. The free energy becomes a func- 
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tional of these density functions, which is claimed to be exact near its min- 
imum. Therefore the variational estimate (w.r.t. density functions), with 
a fixed energy of the system, yields the exact free energy. The string hy- 
pothesis formulation can be regarded as a micro-canonical approach. It is 
supported by many consistency tests. We conclude that within the string 
hypothesis approach the diagonalization is achieved, but the summation is 
cleverly avoided. 

In order to evaluate thermal expectation values of operators, it is better 
to deal with the canonical ensemble. We therefore consider an alternative 
approach based on the Quantum Transfer Matrix (QTM). 4 ' 5 It utilizes an 
exact mapping between a ID quantum system at finite temperatures and 
a 2D classical system. At first sight the formulation may look tautological 
and may seem to be suffering from the need of "summation" . Yet, the main 
claim of the QTM formulation is that this is not the case. As in the the 
string hypothesis approach the "summation" can be avoided. Moreover, 
the QTM makes the evaluation of many quantities of physical relevance 
straightforward. 

This article is organized as follows. In Sec. 2, we present a review on 
the QTM formulation. The results for the bulk quantities will be summa- 
rized in Sec. 2.3. In the rest of Sec. 2, we supplement arguments to justify 
the formula in Sec. 2.3. The non-linear integral equation (NLIE) will be 
explained in Sec. 3 together with an example for the explicit evaluation 
of bulk quantities. The evaluation of the reduced density matrix elements 
(DME) will be discussed in Sec. 4. 

2. The QTM formulation 
2.1. The problem 

Let H be the Hamiltonian of a ID quantum system of size L and V its 
space of states. Our goal is to calculate the thermal expectation value of 
any physical quantity O at temperature T(= l//3) a in the limit L — > oo, 

3 

Here Ej stands for an eigenvalue of H. 

The definition requires both diagonalization and summation. Below we 
shall show how we can avoid the latter within the framework of QTM. 



a The Boltzmann constant ks is set to be unity in this report. 
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2.2. The Baxter-Liischer formula 



To be concrete, we specify a Hamiltonian. As a prototypical integrable 
lattice system we choose the ID spin i XXZ model, 

L L 

H = J E(^i+i + + A (<^i+i + !))= E^+i ( 2 ) 



3=1 



where the a a (a = x, y, z) are the Pauli matrices. The periodic boundary 
conditions (PBCs) imply o£ , j = of. The anisotropy is parameterized as 
A = cos 7. The Hamiltonian acts on "the physical space" Vph ys := ®t=i V? 
where V,- denotes the jth copy of a two-dimensional vector space cie + + 
C2e_. The trace in (1) must be performed over V^hys- By definition the 
"Hamiltonian density" /ijj+i is the jth summand in the first sum in (2). 
It acts non-trivially only on Vj (g> Vj+i- 

The above Hamiltonian is integrable in the following sense. Let R(u, v) 
be the Ugfel?) R matrix, 6 



R(u,v) 



V 



q u _ q -u 

q - q- 1 



depending on the spectral parameters (or rapidities) u, v € C. We define 
Ep s.t. (Eg)^j = 5 ay i5f3,j- Then the matrix elements R^J can be read off 
from 



R(u,v)= Yl R a plMEi®E 

a, 0,7, 5=1,2 



The index 1(2) refers to e + (e_). See fig. 1 for a graphic representation. We 

8 + + + 



CYi - 



-0 



u 
+->- 



A v 



It 

— ■ 



u 

+ H - 



Fig. 1. A graphic representation for R^2(u,v) and some examples 



put arrows, to distinguish the R matrix from other R matrices appearing 
below. The reader should not confuse them with physical variables. 
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By Rj.j + i(u, v) we mean the R matrix acting non-trivially only on the 
tensor product Vj(u) <g> Vj + i(v) of C/ g (s [2 ) modules. We also introduce the 
intertwiner j +1 (u,v) = Pjj + iRjj + i(u,v), where P : Vj(u) ® Vj+i(v) — > 
Vj+\(v) (g) Vj(u). Then, with q = e* 7 , we have the expansion 

Rl 3+1 (u, 0) = 1 + ^j^—u (h jtj+1 + h' j>j+1 ) + 0(u 2 ), 

where h'jj+i := iJsin7(cr| — We introduce the row-to-row (RTR) 

transfer matrix Irtr(w) E End(V^,hys), 

Trtr(-") = tr a R aiL (u,0)R atL -i(u,0) ■ ■ ■ R a ,i(u, 0). (3) 

With the lattice translation operator e lP , shifting the state by one site, we 
obtain the Baxter-Liischer formula 3 

T BTR (u)=e* p (l + —^H + 0(u 2 )). (4) 
v 4 J sin 7 ' 

Note that the h'jj+i terms cancel due to the PBCs. The huge symmetry 
U q (sfa) is at the bottom of the integrability of the Hamiltonian. 



2.3. A summary of results for bulk quantities 

We first present the formula for the free energy per site in the QTM for- 
malism. A supplemental discussion will be given in subsequent sections. 

We introduce the transposed R matrix Rj k (u,v) 7 by (R t )JJ(u,v) = 
R Sa a{v,u). See fig. 2. The QTM does not act on V^hys but on a fictitious 



5 

u ii 



+-► 



U <7 



Fig. 2. A graphic representation for (B*)o 7 d) and some examples 



space Vrrotter = V\{u) ® V2{— u) ® ■ • • Vn-i(u) <E> Vjv(— u). The fictitious 
system size N is often referred to as the Trotter number. The parameter u 
is fixed to be 

4/3Jsin7 4Jsin7 
U ~ jN ~ jNT 
In its most sophisticated version, the QTM is explicitly defined by, 7 
Tqtm(x, u) = tY a R aN {ix, -u)Rl N -i(ix, u) ■ ■ ■ R a 2(ix, ~u)R f al (ix, u). (5) 
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The new parameter x will later play the role of a spectral parameter. The 
factor i is introduced for convenience. 

We are now in a position to write down the formula for the free energy 
per site in the thermodynamic limit, / = — linir J _ ) . 00 lnZir>(P)- 

Theorem 2.1. Let A be the largest eigenvalue of Tq TM (0,u). Then the 
free energy per site is solely given by Aq, 

f = - lim Tin A . (6) 

The limit N — > oo is referred to as the Trotter limit. As was announced 
earlier, eq. (6) expresses / without recourse to any summation. We also note 
that In A itself is already intensive, which may reflect the size dependent 
interaction of the system. 

The quantitative analysis of (6) is most efficiently performed by means 
of the NLIE. Having in mind the examples, from now on we are considering 
only for J — ^,7 — > and u = — jjj, consequently. Let a be the unique 
solution to the NLIE b 

In a(x) = /3e (z + i)-J c {x _ 2 y)2 + A In2l0/)f (7) 

^ x) = h+ { X -i)( X + i) * :=1 + a - 

Here the contour C is a closed narrow loop which encircles all "Bethe 
roots" . We added a Zceman term | ■ ct| to the Hamiltonian so that 
diag(exp(— exp(^)) is inserted in the trace in (5). Then we have the 
following 

Theorem 2.2. The free energy per site can be evaluated in terms of the 
solution to the NLIE. 

" / = ! (1+k) -/ c :RjW" a(l) v- (8) 

Note that the NLIE (7) and the expression for / in (8) are independent 
of N. The extension to arbitrary J, 7 is straightforward. 

Below we shall comment on the derivation of the formula. By presenting 
supplementary arguments, we wish to convince the reader that the above 
formalism, seemingly complicated, is actually necessary and efficient for 
many purposes. Hereafter we set again h = for simplicity. 



b To be precise, there are, in general, several equivalent versions of NLIEs. We present 
one of these below. 
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2.4. The ID quantum partition function as a 2D classical 
partition function 

We define a rotated R matrix R(u,v) by R^](u 7 v) — R]^(v,u) (fig. 3). 
Then we introduce a rotated transfer matrix Trtr(m) S End(Vph ys ) by 

Trtr(u) = tT a R a x(-u,0)R a ,L-i(-u,0) ■ ■ ■ R a , 1 (—u, 0) . 

Analogous to (4) Tr TR (m) = e" lP (l + uR + 0(u 2 )). We thus obtain an 
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Rp](u,v) 



v—u 
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v: — a 
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Fig. 3. A graphic representation for R^g(u,v) and some examples 



important identity, 

Zi O {0) = tr yphy8 e"^ = Jim tr yphy3 (T doublc ( U ))f | u ^_ » (9) 

where Td ou bie(M) := 7rtr(u)2rtr(m). The rhs of (9) can be interpreted as 
a partition function of a 2D classical system defined on N x L sites (fig. 4), 

Zir>(/3) = lim Z 2 Dciassica.i(N, L,u — --rz). 

N—tcc iv 

This equivalence lies in the heart of the QTM formalism. The expression 
(9) itself, however, is of no direct use for the actual evaluation of physical 
quantities for the following reason. Let the eigenvalue spectrum of Trtr( u ) 
be Xq(x) > Xi(x) > X 2 (x) > ■ ■ ■ . We introduced x = i~ 1 (u+l) for technical 
reasons. It is easy to see that Tr TR (u) has the same spectrum. Thus, 

tr w ( Wu)) * = (AoW) N (l + (^g) N + N + • • • ) ■ (10) 

The eigenvalue Xj (x) is characterized by its zeros ±0 a (a = 1, 2, . . . ) on the 
real axis (holes). We know numerically that for low excitations, 9 a ~ hiL 
and also that Xj(x) is analytic and nonzero in the strip |S&ria;| < 1 except 
at ±9 a . Let us introduce an analytic and nonzero function near the real 
axis, Xp by Xj(x) = ]la tn f ( x ~ ®a)thj(x + 9 a )X^(x). It approximately 
satishes the inversion relation for L ^> 1, 

\\{x-i)\]{x + i) = 4>(x), (11) 
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Fig. 4. Fictitious two-dimensional system 



where 4>{x) is a known function common to any j. Thus, we simply have 



Xj(x) 




X (x) 





9 a )th-(x + e a ) 



For very low excitations, we take a single pair of holes, substitute 8 a ~ 
~ In and take the large L limit. Then we arrive at the estimate [u ~ 0) 



\j{x) 



X (x) 



thus 



Aj (x) 



\ (x) 



N 



|u| 



(12) 



For a usual 2D classical system we can consider an infinitely long cylinder 
and take ^ S> 1. We thus have to take into account only the first term 
on the rhs in (10). By contrast, the spectral parameter depends on the 
fictitious systen 
T ^ 0, we have 



fictitious system size u = — in the present case. Therefore, as long as 



\j{x) 



X (x) 



N 



O(l) forL>l. 



Fig. 5 presents numerical evidence for the above argument. The left figure 
shows the histogram of the distribution of | A^- / Ao | for q — 1,L — 10, u = 
—0.01 in the sector with vanishing magnetization. One clearly sees that 
the maximum of the distribution lies near (A^/Aol ~ 1. The right figure 
magnifies the region near | A^- / Ao | ~ 1- The maximum is located around 
(Aj/Aol ~ 0.96. We believe that, with increasing L, the peak moves towards 
\Xj/\o\ ~ 1. These findings are consistent with (12). Hence, we conclude 
that infinitely many terms of the sum in the rhs of (10) contribute non- 
trivially, and eq. (9) is of no practical use. 
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oRow spectra; L=10u=-0.01 



Fig. 5. The distribution of eigenvalues. The horizontal axis is the absolute value of 
the eigenvalues normalized by the largest one. The left figures ranges over [0,1] in the 
horizontal direction, and the right one is zoomed into the range [0.93,1], 



2.5. Commuting QTM 

A crucial observation was made in reference. 4 We start from the same two- 
dimensional classical model in fig. 4. We consider, however, the transfer 
matrix propagating in horizontal direction, that is, Tq TM (u). Equivalently, 
one can rotate the system by 90°. Then we define a transfer matrix prop- 
agating in vertical direction, Tqtm(u) (see fig. 2.5). The latter is more 
convenient for our formulation. 

The partition function is then given by, 

Z m {P) = Jim try Trottcr (T QTM (u)) L (13) 

TV— J-OO U — —-TT 




L 

Fig. 6. The graphical definition of Tqtm- 



Let the eigenvalue spectrum of Xqtm(w) be A (u) > Ai(tt) > A 2 (u) > 
■ ■ ■ . Then we have an expansion similar to (10) 

trv Trotte {Tqtm{u)) L = (A («)) i (l + (^) L + (j^) L + •••)■ (14) 
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Our physical interest is in the free energy per site / in the thermody- 
namic limit L — > oo. 



f = — — lim lim 

jj L— s-co JV^co 



1 



■In 



{lnAo(w) 



\A a (u)J J J 



L V \A (u) 
Proposition 2.1. The two limits in (15) are exchangeable 



(15) 



We supplement an argument which claims that the second term in the sec- 
ond line in (15) is negligible for L — > oo. The previous argument, using the 
inversion relation (11) can not be applied directly as the spectral parameter 
u is already fixed as — ^ in the present problem. 

We introduce a slight generalization, a commuting QTM Tqtm(x,u), 
by assigning the parameter ix in "horizontal" direction. 10 The substitution 
x = recovers the previous results. The precise definition is shown in (5). 
Hereafter we drop the u dependence as it is always — jj . Let 7qtm (x) be the 
corresponding monodromy matrix. Then it is easy to see that monodromy 
matrices are intertwined by the same R matrix as in the RTR case, 

R(x,x')Tqtm(x) <8>7qtm(z') = 7qtm(2;') ®Tqtm{x)R(x,x'). (16) 

This immediately proves the commutativity of Tqtm(x) with different x's. 

The most important consequence of introducing x is that we have the 
inversion relation in this new "coordinate" , 

A i j (x-i)A* j (x + i) = i>(x,u) (17) 

where we set again Aj(x) = J| a th j(x — 9 a ) th j(x + 9 a )A*-{x). Note that 
A!- (a;) also depends on the "old" spectral parameter u, which is set to be 
— jj on both sides. The known function ip is again independent of j. The 
analysis of the Bethe ansatz equation associated to the QTM implies that 
a ~ - In 4^- for large /3. Then, proceeding as before, we obtain, 



Aj(x) 



A {x) 



- ch %x 



thus 



Aj(x) 



A (x) 



~ e 



- ch f x 



The diagonalization for fixed N clearly shows the gap between the eigen- 
values, which is consistent with the above argument. Thus, at any finite 
temperature, the second term in (15) is negligible for L — > oo. We then 
conclude that the formula (6) is valid. 

Although we made use of the integrability of the model in the above 
argument, the validity of the formula is actually independent of it. See the 
proof in reference. 4 
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3. Diagonalization and NLIE 
3.1. Bethe roots 

Thanks to (16), one can apply the machinery of the quantum inverse scat- 
tering method, devised originally for the diagonalization of Trtr, to the 
diagonalization of Tqtm- We skip the derivation and present only results 
relevant for our subsequent discussion . We fix N for a while. Then the 
eigenvalue of Tqtm is given by 

a(x) := <j)+{x + 2i)cj)-(x) d(x) := <p-(x — 2i)(p+(x) 

m , . 

Q(x):=Y[(x- Xj ) ^ ±(a;):= (™)-. 
j'=i 

The different sets of Bethe roots {xj } correspond to the different eigenval- 
ues. They satisfy the Bethe ansatz equation (BAE), 



a(xj) Q(xj + 2i) 
d(xj) Q(xj — 2i) 



1 < j < m. (19) 



For the largest eigenvalue the number of roots m equals -y. 

To evaluate / via (6) we need the largest AW for N — > oo. This means 
we must deal with infinitely many roots in the limit, which resembles the 
situation encountered in the evaluation of the free energy in the thermody- 
namic limit of a classical 2D model by means of the RTR transfer matrix. 
Still, we would like to comment on the qualitative difference in the root dis- 
tribution between such "standard" case and the problem under discussion. 

Fig. 7 shows the distribution of the positive half of BAE roots for the 
largest eigenvalue of Trtr (left) and Tqtm (right) for various system sizes. 
The distribution of RTR roots behaves smoothly for large system size. The 
limiting shape of the distribution (the root density function) is a smooth 
function satisfying a linear integral equation. For Tqtm, on the other hand, 
a few large roots remain isolated at almost the same positions as N in- 
creases, while close to the origin more and more Bethe roots cluster. 

Let us describe this in detail. Using the NLIE technique, we can derive 
an approximate BAE equation (see the discussion after (25)), 

AT^. thf^.-^ x 



C A technical remark: the vacuum is conveniently chosen (+, —,+,—,•••). 
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KLilk-iiiil R,m m R,r\ BAH Run], L=K. I 6,32.64. 1 23 



Rational QTM BAB RihiIs lx'Li=i>.4 N=3.16.32,64. 123 



Fig. 7. The positive half of BAE roots for RTR (left) and for QTM (right) with system 
size, from 8 (bottom) to 128 (top). 



The hole 9 a corresponds to the branch cut integer I max = x — §; an d this 
implies 9 a ~ \ ln/3 for j3 ^> 1, which was used in the last subsection. 

— A~ 

c 3 ~ N X 3 



Near the origin, we set Xj = jrXj, and obtain the approximate distri- 



bution of Xj as an algebraic function, 



p{x)= lim Ij - 



jv^oo N(x j+1 - Xj) 2n(x 2 + 1) 

This differs from the usual root density function which decays exponentially 
as \x\ — > oo. In the original variable, if we take the Trotter limit naively, 

p{x) ~\% M^T u*) ~\ 5{x) - 

Namely the distribution of the BAE roots for Tqtm is singular in the Trotter 
limit. We thus conclude that the usual root density method may not be 
applicable, and we have to devise a different tool. 

Let us stress again that the cancellation of of the order of N many terms 
in lnA^' is a unique property of the QTM. In the RTR case O(L) many 
terms can survive, and we obtain intensive quantities (e.g., the free energy 
per site) only after dividing by L. On the other hand, In A( w ) is already an 
intensive quantity, as remarked after Theorem 2.1. The cancellation is thus 
vital. The 0(e N ) terms, a(x),d(x) must be canceled by the denominator 
Q(x), resulting in 0(1) quantities. According to this point of view, (18) is 
(again) not practical, as the ratios of 0(e N ) terms are still present. Thus, 
we understand (18) still as a starting point, not as the goal. 

One must intrinsically deal with a finite size system with coupling con- 
stant u depending on the Trotter size, and then take the Trotter limit. Such 
an attempt was executed first numerically 8 by extrapolation in N . The ana- 
lytic low temperature expansion was performed 9 based on the Wiener- Hopf 
method. Below we shall present the most sophisticated approach which uti- 
lizes the commuting QTM in a most efficient manner. 7 
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3.2. Non-linear Integral Equation (NLIE) 

The introduction of the new spectral parameter x plays a fundamental 
role. Instead of dealing with the BAE roots directly, we make use of the 
analyticity of specially chosen auxiliary functions in the complex x plane. 

There are various approaches. One of them is to introduce the fusion 
hierarchy of the QTM, which contains the original Tqtm as T\. In place of 
the BAE one uses the functional relations among fusion transfer matrices, 

T m (x - i)T m {x + %) = i> m (x) + T m -i(x)T m+ i(x) . (21) 

Here tpi is nothing but ip in the inversion relation (17), where the small 
term TqT 2 /ipi was neglected. As {T m } constitutes a commuting family, 
the same relation holds among the eigenvalues. We thus use the same 
symbol T m for the eigenvalue. After a change of variables, y m (%) — 
T m _i(x)T m+ i(x)/il) m (x), one can transform the algebraic equations into 
integral equations under certain assumptions on the analyticity of y m . The 
resultant equations coincide with the Thermodynamic Bethe Ansatz equa- 
tions. 10,11 The string hypothesis is thus replaced by an assumption on the 
analyticity of y m . The coupled set of equations may fix the values of T\(x). 
Then A = Ti(0) yields the free energy per site /. A technical problem in 
this approach is that we must deal with an infinite number of y m functions, 
which requires a truncation of the equations in an approximate manner. 

Below we shall discuss another approach originally devised in the con- 
text of the evaluation of finite size corrections. 12 We define the auxiliary 
function <Xn(x) by the ratio of the two terms in \( N \x) (18), 

. . d(x) Q(x + 2i) . . , . 

a ^-^)W^ *M*) = i + M«). 

The suffix N is introduced to recall that we are fixing N finite here. The 
BAE (19) is equivalent to the condition 

aN{xj) = — 1 or In dN(xj) = (2/j + l)ni. (22) 

We also note that lim^i^^ cin(x) = 1 by construction. 

We then adopt the following assumptions for the analytic properties 
of 2ljv(#) corresponding to the largest eigenvalue. They are supported by 
numerical calculations. 

(1) There are ^ simple zeros of %In{ x ) on the real axis. They coincide with 
the BAE roots. There are additional zeros, sufficiently far away from 
the real axis, so that C does not include them inside. 
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(2) The only pole of QLn{x) in e [—1, 1] is located at x = iu and is of 
order ^. 

Once these assumptions are granted, one immediately derives the fol- 
lowing NLIE, 

/ x , d>-(x + 2i)6+(x) f 2 , „, . .dy 

4> + (x + 2l)(t)-{x) Jc(X-y) +4 7T 

The largest eigenvalue A can be similarly represented by 

>»A""W=>»M* + 2.)*-(, - 2i)) + / c I _|*<?L_|, (24, 

Note that only the driving term in (23) depends on N. We can thus take 
the Trotter limit easily, with a := limAr->.cx> ciat, and obtain the NLIE in (7) 
(for h = 0). To evaluate the free energy one has to first set x = 0, then take 
the Trotter limit. Or otherwise one meets a spurious divergence. Then we 
obtain the expression for the free energy in Theorem 2.2. 

One still needs to make an effort to achieve a high numerical accu- 
racy, especially at very low temperatures. The introduction of another 
pair of auxiliary functions solves this problem. We define ajv,2tjv by 
a N (x) = (ojv(ar)) -1 , %n(x) = 1 + a N (x). 

Numerically one finds that |a/v| ;s 1 for Qx ^ 0. Thus, we use a at, 21 at in 
the upper half plane and Oat, 21 at in the lower half plane. It is straightforward 
to rewrite (23) in the coupled form, 



lna N (x)=Df\x)+J c F\x~y) lnS^^-jf F\x~y) h& N (y)^ (25) 

\na N (x)=D ( _ N) (x)+f F{x-y)m N ( y )p- [ F{x-y)\m N {y)^- 
Jc- 2n Jc+ 2n 

± 2 v thf(a;-m) ; ' v ; J^l + e 2 ^ ' 



where C+(C_) is a straight contour slightly above (below) the real axis. In 
the first (second) equation we understand that x G C+(C_). Note that the 
convolution terms bring only minor contributions as they are defined on 
those contours where the auxiliary functions are small. Therefore the main 
contributions come from the known functions. This enables us to perform 
numerics with high accuracy. We can drop the convolution terms for the 
lowest order approximation. Thanks to eq. (22) this leads to eq. (20). 
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Similarly, for the largest eigenvalue we have 

lnAW(j)=eW(i)+/ K + (x - x)\nSl N (x')—+ [ K-(x - x)ln$l N (x')— , 
Jc+ 2n JC- 2lT 

K ± (x)=K(x±i), K(x)= J 

£ W = ln0 + (x + 2iU(x - 2i) - ? [e-^ kx ^dk. 

2 J kchk 

We obtain the NLIE and the eigenvalue in the Trotter limit by replacing 
Oat — > a etc. and 

± 2shfx 2 V J l + e 2 l fe l ' 

For the actual calculation, it is even better to deal with b(x) := a(x + i) 
and b(x) :— a(x — i) so that the singularities of ln(l + b), ln(l + b) are away 
from the integration contours. We omit, however, the details. 

As a concrete example for the evaluation of bulk quantities we plot 
the susceptibility, \ — d^f, in fig. 8 (left). Note that at low temperatures 
the physical result in the Trotter limit (solid line) deviates from its finite 
Trotter number approximation. The above approach has been successfully 




Fig. 8. Left: the susceptibility of the s = i XXX model, in the Trotter limit (solid line) 
and for fixed Trotter number (crosses:^ = 128, squares: N = 1024). Right: a plot of £T 
against temperature for the XXZ model. 



applied to many models of physical relevance. 13-19 

The correlation length £ characterizes the decay of correlation functions 
at large distance, e.g., 

(a+a-)~e-V L (26) 
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It is evaluated from the ratio of the largest and the second largest eigen- 
values of the QTM. 9 ' 20 ~ 22 For the second largest eigenvalue state, our as- 
sumption (1) on 2t is no longer valid: a pair of holes 8 a lies on the real 
axis, and they are zeros of 21 other than Bethe roots. Nevertheless, a small 
modification leads to a set of equations that fixes the second largest eigen- 
value. The resultant NLIE has a form similar to (25) containing, however, 
additional inhomogeneous terms. Fig. 8 (right) shows a plot of £T against 
temperature for the XXZ model with q = c l ~ for p = 3,4, 5. 21 

When ft 7^ 0, we have to replace a(x),d(x) in (18) by e~ tjh / 2 a(x), 
c' 3h / 2 d{x). Then we add /3ft to the rhs of (23) ((24)). Also D { ± } 

must be replaced by ± 

Before closing this section, we would like to mention another formulation 
of thermodynamics also based on the QTM. 23 It is described by a NLIE 
for Aj directly and allows one to efficiently calculate high temperature 
expansions. The good numerical accuracy in the low temperature region 
is, however, hard to achieve. Moreover, we point out that the equation is 
the same for any eigenvalue. Thus, one should know a priori good initial 
values in order to select the convergence to the desired eigenvalue. 

4. DME (density matrix elements) at finite temperatures 

The deep understanding of a model requires ample knowledge of its cor- 
relation functions. We would therefore like to go beyond their asymptotic 
characterization by the correlation length £ (26). 

The evaluation of correlation functions has been defying many chal- 
lenges in the past. Considerable progress was made only recently for the 
T = correlations, based on vertex operators, 24 on the qKZ equation 25 
and on QISM. 26 The third approach is the most relevant for our purpose. 
For T = it first requires the solution of the "inverse problem", that is, 
one has to represent the spin operators in terms of the QISM operators 
A{u),B{u),C{u),D{u). Then, by algebraic manipulations, one obtains the 
correlation functions as combinatorial sums of expectation values of QISM 
operators, which are finally converted into (multiple) integrals. 

At first glance, the case T > seems far more difficult, as one expects 
that a summation of the contributions from all excited states is necessary. 
We argue here that, as above, the QTM helps us to avoid this summation 
and that, moreover, one docs not have to solve the "inverse problem" within 
the QTM framework. The combinatorics, on the other hand, can be done 
in parallel to T = 0, because the QISM algebra is the same in both cases. 

Let us explain why we can bypass the inverse problem in the QTM 
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formalism. This can be most quickly done in a graphical manner. To be 
specific, we need to evaluate DME, 

tr Vh e-WEZ 1 ■ ■ ■ E%™ 
U l -p m - ■ ■ ■ &p m > ~ tr Vphy .e-^« ' 

Using the logic of section 2, we can represent e _/3w by a "2D partition 
function". Therefore -D^ 1 ...^" can be represented by a modified 2D parti- 
tion function: Start from the N x L classical system (fig. 4) with periodic 
boundaries in both directions. Then cut n successive vertical bonds at the 
bottom row, and fix the variables at both sides of the cut. As we are adopt- 
ing PBCs in the vertical direction this is equivalent to fixing the configu- 
ration of n successive bonds at the top and at the bottom. See fig. 9 (left). 




As previously, we rotate the lattice by 90°. See fig. 9 (right). Obvi- 
ously we can write -D^ 1 ...^™ in terms of elements of the monodromy matrix 
7qtm(^)- By introducing independent spectral parameters £j we obtain 

(r>Y U '"' an (t \ = ^01(^)^(6) ■ ■ ■ (T Q TM)g:(£n)l<&0) 

V 1 h,-A 1 ' ($o|rQT M (6)---r Q TM(^)|*o> ' 

Here <&o denotes the largest eigenvalue state of the QTM, which is given by 
acting with B operators on the vacuum. As any (7qtm)^ is represented 
by a QISM operator, we reach an expression for DME purely in terms of 
QISM operators without solving the inverse problem. 

At the same time, the problem for T > is not so simple in view of the 
analyticity. We consider D++ as a concrete example. After employing the 
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standard QISM algebra, one obtains 

Dttitub) x 21(6)21(6) = (E j>fe 



(gfc-|2)faj-|i-2t) 



2+2i 



6,2 ^3 (x 3 -£ 2 +2i) W J:t £ 2jl (z.,-£i+2t) "O,: 

Here a;j denotes a BAE root and 21 is the auxiliary function. We introduced 
Wj t k in order to deal with the ratio of inner products of wave functions. 
Wj.k is characterized by a simple algebraic relation. Note that the above 
algebraic expression is formally identical for T = and T > 0: one only has 
to replace Xj and Wj.k for T = by those for T > 0. 

In the case T = there are several simplifications. First, the auxiliary 
function is by construction trivial, 01= 1. Second, we can introduce the root 
density function in the thermodynamic limit. Then the algebraic relation 
for w jtk -> g{xj,ik) is solved with the explicit result g(x, £) = 4ch ^ ( ^._^ + - ) . 

g{x,&) g{x',&) 



l+2i 



-x k -2i) 
>j-gi) 



Wj,2 Wfe,2 
2 + 1 ) • 



2,) 



^ll(a,6) = (/^ ( tffc"l 1 2, ) 



;g{x,&) + 1 



6,2 J M ""(x-c 2 +2i)y^7Si; £ 21 j "^( x _^ 1+2 i): 

Third, we can freely move the integration contours. Every time it passes 
the singularity of g(x), it brings extra contributions and they cancel the 
"tails" (the 2nd to the 4th terms above). We finally obtain 

g(x+i,£i) g(x'+i,£i) 
g{x+i,&) g(x'+i,&) 



D ++ 



(6,6)= 



oo ^^-^ +0(x " gi r ) 

-oo b,i{x-x ' -2i) 



(27) 



Without such a compact expression, it is hard to proceed further. 

On the other hand, 21 is quite non-trivial for T > 0. As noted previously, 
we can not resort to the root density function. The explicit form of Wj t k in 
the Trotter limit is thus unknown. The most significant difference is that 
the integration contour is already fixed for T > 0. Thus, we cannot apply 
the above trick to swallow tails into the ground state. 

Nevertheless, with an appropriate choice of a further auxiliary func- 
tion G(x,£), it was shown that a compact multiple integral representation, 
similar to (27) is also possible for T > 0, 27 ' 28 



dx' (jc-6 -2i)(x'-&) 



G(*,£i)G(x',£i) 
G(x,6)G(x',6) 



Ic 2l(x') 4ir 2 £ h2 (x - x' -2i) 

The formula for any other DME is similarly known. 

It is a big progress to obtain the multiple integral representation for 
DMEs. The representation is, however, not yet optimal. Although one can 
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use it for the numerical analysis at sufficiently high temperatures, it suffers 
from numerical inaccuracy at low temperatures. 29,30 We thus would like to 
reduce it to (a sum of) products of single integrals. 

The factorization of DME at T — has been performed by brute force, 
with the extensive use of the shift of contour technique. 31,32 Based on stud- 
ies of the qKZ equation, a hidden Grassmannian structure behind DME 
has been found. 33 ft naturally explains the factorization of the multiple 
integral formula through the nilpotency of operators. The explicit form of 
DME consists of two pieces, the algebraic part, evaluated from a matrix el- 
ement of g-oscillators, and the transcendental part, related to the spinon-5 
matrix. 

On the other hand, we still do not have a finite temperature analogue 
of the qKZ equation. We are nevertheless able to factorize the multiple 
integrals for small segments. 34 The explicit results also consist of two parts. 
Surprisingly the algebraic part remains identical to the T = case, while 
the transcendental part can be interpreted as a proper finite temperature 
analogue to the spinon-S 1 matrix. This finally enables us to perform an 
accurate numerical analysis of the correlation functions. 35 We show plots 
of (crf<j|) for A = -4= with various magnetic fields in fig. fO. The results 
from a brute force calculation are also plotted, which supports the validity 
of our formula. Such high accuracy calculation can clarify the quantitative 




nature of interesting phenomena such as the quantum-classical crossover. 36 
Recently a proof of the existence of factorization of the DMEs for T > 
was obtained, again by using the Grassmannian structure. 3 See also the 
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further development 38 in this direction based on NLIEs. 
5. Summary and discussion 

We presented a brief review on the recent progress with the exact thermo- 
dynamics of ID quantum systems. The QTM is found to be an efficient 
tool, and it offers a framework to evaluate quantities of physical interest, 
including DME. The NLIE combines into the framework nicely, yielding 
high accuracy numerical results. 

The factorization of the multiple integral formula at T > is yet to be 
further explored. It seems e.g. quite plausible that the qKZ equation could 
be extended to finite temperatures. This might be an important next step. 

There are certainly many interesting questions left open. For example, 
can we have the QTM formulation starting from a continuum system? What 
is the generalization of the multiple integral formula to models with higher 
spin? The study of such questions is underway. 
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